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Abstract 

. . .Acoustic neutrino detection is a promising approach for large-scale ultra-high energy neutrino detectors in water In this article, a 
f^ Monte Carlo simulation chain for acoustic neutrino detection devices in water is presented. It is designed within the SeaTray/IceTray 
^«. software framework. Its modular architecture is highly flexible and makes it easy to adapt to diff'erent environmental conditions, 
*vj detector geometries, and hardware. The simulation chain covers the generation of the acoustic pulse produced by a neutrino inter- 
j^ action and the propagation to the sensors within the detector In this phase of the development, ambient and transient noise models 
O I for the Mediterranean Sea and simulations of the data acquisition hardware, similar to the one used in ANTARES/AMADEUS, are 
.^ implemented. A pre-selection scheme for neutrino-like signals based on matched filtering is employed, as it can be used for on-line 
filtering. To simulate the whole processing chain for experimental data, signal classification and acoustic source reconstruction 
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algorithms are integrated. In this contribution, an overview of the design and capabilities of the simulation chain will be given, and 
1 some applications and preliminary studies will be presented. 
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Q 1. Introduction 

Tj^ Acoustic neutrino detection uses the effect that ultra-high en- 
Cu^ ergy (UHE) neutrinos can produce a detectable acoustic pulse 
according to the thermo-acoustic model yj]. This model de- 
scribes the generation of an acoustic pulse due to the local heat- 
ing of the medium by a hadronic shower, which is caused by 
an UHE neutrino interaction. The fast deposition and slow dis- 
sipation of the energy by the cascade in the medium leads to 
a bipolar pulse (BIP) and due to the cylindrical geometry of 
the shower the wave propagates through the medium in a disk- 
like shape perpendicular to the main axis of the cascade. Given 
fT^ the expected low flux of neutrinos with energies in excess of 
T-H 1 00 PeV, a potential acoustic neutrino telescope must have large 
'>^ dimensions of presumably > lOOkm^. The attenuation length 
• 'H of sound in water is of the order of 1 km for the peak spectral 
rN density of around 10 kHz, allowing for a less dense instrumen- 
tation. In this article, the efforts taken to develop a complete 
" simulation chain to reproduce the acoustic pulse generation, the 
detector properties, and the deep-sea acoustic environment are 
reported. 



2. Simulation Chain 

In this section, the design of the simulation chain and the ca- 
pabilities of the different modules will be presented. The simu- 
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lation chain consists of the following stages, which build up on 
each other to create a simulated event: 

• An interaction vertex is located at a random position in 
a given volume around the detector and the energy and 
direction of the shower are set randomly. 

• The formation of the shower and the resulting acoustic sig- 
nal, which is generated by an UHE neutrino interaction, 
are simulated. 

• The acoustic environment of the deep sea is reproduced 
including both the ambient and transient noise conditions. 

• The data acquisition (DAQ) hardware is simulated includ- 
ing the system response and inherent noise of the sensors 
and read-out electronics. 

• A pre-selection scheme is applied. 

The simulation chain was designed within the SeaTray/IceTray 
software framework |J2l |3|]. Its modular architecture is highly 
flexible and makes it easy to adapt it to different environmental 
conditions, data acquisition hardware, and detector geometries. 
The generation of the Monte Carlo (MC) shower is the first step 
described here. The MC shower is produced from a parame- 
terization, which is based on work by the ACoRNE collabora- 
tion [4, 5]. This parameterization describes the distribution of 
the deposited energy in the surrounding medium, in this case 
water. The longitudinal and radial energy density distribution 
of a 10" GeV shower is shown in Fig.[T] From this energy dis- 
tribution, a point distribution is produced with a point density 
proportional to the energy density distribution. The represen- 
tation of the shower is shown in Fig. |2j it is about 20 m long 
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Figure 1 : The longitudinal and radial distribution of the deposited energy of the 
shower in water. The total shower energy is 10 " GeV. 
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Figure 2: MC point distribution generated proportional to the deposited en- 
ergy distribution of a lO" GeV hadronic shower as in Fig. [T] Bin sizes are 
0.001 m in x and 0.005 m in z, where the points have been projected from a 
three-dimensional distribution upon the .r^-plane. 



and has a diameter of about 10 cm. After the cascade has been 
simulated, the acoustic pulse and its propagation to the sensors 
within the detector are calculated. The deposited energy of the 
shower produces a local heating of the medium. With respect to 
hydrodynamical time scales, the energy deposition at time to is 
instantaneous and the dissipation of the energy is slow in com- 
parison. The energy deposition e(r, t) can be factorized into a 
spatial and temporal part using the Heaviside function. 



e(r,f) = e(r)0(f-fo) ^ -e(r,t) ^e(r)6(t-to) (1) 
ot 

Assuming a total energy deposition E with cylindrical symme- 
try, the spatial part e(r) can be expressed for the longitudinal 
and radial positions z and r in the shower as: 



e(r) 
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Inserting Eq.[T]and Eq.|2]into the expression for the pressure p, 
following ||6D: 
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where a is the thermal expansion coefficient, Cp is the specific 
heat capacity at constant pressure and Cs the speed of sound in 
the medium. Eq. [3]can be further reduced to: 



p{r, t) 



Ea 

4^ 



r dh' , 

e(r' 



dt Cs 



(4) 



where /? = |r - r'| is the distance between the shower maximum 
and the sensor. A velocity potential can be defined: 
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As described before, the distribution of points within the MC 
shower is proportional to the energy density distribution of the 
shower. So the pressure at a sensor in a distance R from the 
shower can be numerically calculated. The signal propagation 
time from each point within the MC shower to a sensor in the 
detector is calculated and entered into a histogram with a bin- 
width according to the sampling rate (here 1 jjs) and a size big 
enough to hold the distribution (here 2'^ bins). After normalis- 
ing each bin with the number of points in the shower, scaling it 
with the contant term (see Eq.|5]l, and dividing by the mean dis- 
tance to the shower, this results in the velocity potential E„..(t) 
as shown in Fig. [3] The Fourier Transform of Eq. ^including 
Eq.|5]can be written as: 



p(u) 
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J di 



E,,y,(t)e'"^'dt = iuE„.,(aj), 
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taking into account the basic property of the Fourier Transform 
that the derivative in the time domain is the same as multiplying 
by icij in the frequency domain. Exyz{cL>) is derived form the his- 
togram Exyzit) using the Fast Fourier Transform (EFT). The fre- 
quency dependent sound attenuation in sea water is multiplied 
with the pressure signal in the frequency domain. The attenua- 
tion is based on a model by Ainslie and McColm [7] extended 
to its complex representation. This procedure leads to the char- 
acteristic acoustic signal - a bipolar pulse - at the sensor as 
shown in Fig.lD 

The background for acoustic neutrino detection in the deep 
sea consists of two different types of noise: transient and ambi- 
ent noise. Transient noise signals are short, but of high ampli- 
tude, and can mimic bipolar pulses from neutrino interactions. 
In the simulation, four types of transient signals are imple- 
mented so far: bipolar and multipolar pulses, sinusoidal signals, 
and signals with brown noise frequency characteristics. Sources 
of these four types can be anthropogenic sources like shipping 
traffic or marine mammals. The ambient noise is mainly caused 
by agitation of the sea surface |8], i.e. by wind, breaking waves, 
spray, and cavitations. Thus it is correlated to the weather con- 
ditions, mainly to the wind speed. The model used for the sim- 
ulation of the ambient noise is based on the so-called Knud- 
sen spectra H, which are adapted to the deep sea by applying 
attenuation effects. The power spectrum density (PSD) of the 
ambient noise is shown in Fig. |5l the scatter plot presents the 
PSD of ambient noise for different levels of the wind speed 
and shipping traffic. The wind speed distribution included in 
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Figure 3: The velocity potential Sn-zCf) as a function of the signal propagation 
time, which is normalized to the mean propagation time. The total energy de- 
position of the shower is lO" GeV at a distance of 1000m from the sensor 
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Figure 4: Simulated acoustic pressure signal for a lO" GeV shower as recorded 
at the position of the sensor The distances between shower and sensor is 
1000 m. The bipolar shape of waveform with an asymmetric tail is recognis- 
able. 



the simulation corresponds to measurements at several weather 
stations near the coast of Marseilles, France. The mean noise 
level {(Tnoise) IS about 25 mPa for the frequency range from 1- 
-100 kHz and for 95 % of time the noise level is smaller than 
2{cr„oise} 1 101 ]. This is reproduced by the simulation as shown in 
Fig.lSl 
The simulation of the DAQ hardware comprises two parts: 
the simulation of the sensor and the read-out electronics re- 
sponse. The design of the DAQ hardware is inspired by the 
AIVIADEUS lull project. This includes acoustic sensors using 
the piezo-electric effect (hydrophones) and read-out electron- 
ics to amplify and digitise the signal. The inherent noise and 
the system transfer function, for both the sensors and the elec- 
tronics, have been measured in the laboratory. Sensors normally 
show a directional dependency of their sensitivity, therefor sig- 
nal and ambient noise have to be treated separately. In this case, 
the incident direction of the noise is the sea surface above the 
detector Signal and ambient noise are then superimposed. The 
inherent noise of the sensor is added. For the read-out elec- 
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Figure 5: Power spectrum density (PSD) of the ambient noise as produced by 
the model described in the text. This scatter plot shows the PSD for the different 
levels of the weather conditions and shipping activities. 
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Figure 6: Frequency of occurrence distribution of the ambient noise level as 
function of the ratio between the noise level <T„i,ise of the sample and the mean 
noise level {(Tnoise)- This simulation is based on an ambient noise model as 
described in the text. The left-side y-axis shows the number of occurrences and 
the right-side y-axis the probability that a noise level occurs. 



Ironies, the resulting waveform is convoluted with the system 
transfer function and the inherent noise is added. The output is 
directed to a simulation of an on-line filter system 1.12.] . which, 
for real data, is used to reduce the amount of data to store and 
to pre-select for further off-line analysis. The filter is based on 
a matched filter, which uses a pre-defined bipolar pulse as ref- 
erence to select signals with bipolar shape. In addition, a coin- 
cidence test between the sensors is performed. 

3. Applications 

The simulation chain is used to study signal classifica- 
tion [13J and reconstruction algorithms 1 14]. For the presented 
studies, a detector geometry similar to the configuration of 
AMADEUS fll] is used, assuming six clusters of six sensors 
each at fixed positions and with fixed orientation. The precise 
reconstruction of the arrival time of the signal is crucial for 
the direction and position reconstruction of the acoustic source. 
The arrival time is determined by performing up-sampling of 
the filtered waveform sample and cross-correlation with a pre- 
defined bipolar pulse. This procedure achieves a precision of 
about Ifis. Due to the narrow opening angle of the acous- 
tic emission of a neutrino interaction, local clusters of sensors 
could be preferred in a detector design. Such clusters, consist- 



Deviation form true zenitti 
Deviation from true azimutfi 



Statistics: 
Zenith: 

IVIean = 2.5-10' 

a = 0.70 
Azimuth: 

IVIean = 0.3-10"' 

= 0.69 



120000 

100000 

80000 
w 

■| 60000 

c 

HI 

40000 

20000 



-2-1012 

Angular Resolution ["] 

Figure 7: Angular resolution of the direction reconstruction algorithm shown 
for the zenith and azimuth angle. The mean of the distribution is around zero 
and the sigma is about 0.7 ° for both zenith and azimuth angle. The acoustic 
sources were generated in cube of 5 X 5 X 2.5 km^ around the detector centre. 



Preliminary 




Deviation from true position — 

Statistics; 
Mean = 30.7 
o = 25.5 




20 40 60 

Positioning Resolution [m] 



too 



Figure 8 : Resolution of the position reconstruction of the acoustic source shown 
as the distance between the reconstructed position and the true one. The reso- 
lution peaks at about 5 m, the distribution is broad with a mean of about 30 m 
and a sigma of about 25 m for the given range. The acoustic sources for this 
analysis were generated in cube of 5 X 5 X 2.5 km^around the detector centre. 



ing of several sensors arranged with interspaces of a few me- 
ters, have also advantages for the coincidence test used by the 
on-line filter and for the reconstruction of acoustic source posi- 
tion. The direction reconstruction is based on a least square fit 
of the measured arrival times at a given sensor cluster: 



'<^(t, 



measured^ ^expectedi\"-> y)) /' 



(7) 



where / e I..N (N sensors of the cluster), t the arrival time 
and and (f> the zenith and azimuth angle, respectively. The 
acoustic sources for this analysis were generated in cube of 
5 X 5 X 2.5 km-^ around the detector centre. The angular reso- 
lution reached with this algorithm is centred around zero with a 
sigma of about 0.7 ° for both zenith and azimuth angle as shown 
in Fig. |7] This is consistent with the resolution of the recon- 
struction of the arrival times. The position reconstruction of the 
acoustic source is obtained from triangulation of the previous 
determined direction of the incoming signal. If the directions 
were reconstructed for at least two of the sensor clusters, the 
triangulation is performed by minimising the distance between 
the rays starting at the sensor clusters and pointing into the re- 
constructed direction. The distribution of the distance between 
the reconstructed position and the true one peaks at about 5 m, 
but is also broad with a mean of about 30 m and a sigma of 25 m 
(cf. Fig. [8]l. This is in agreement with the resolution of the di- 
rection reconstruction. 

As a further application, signal classification will now be pre- 
sented. The classification system stems from machine learning 
algorithms trained and tested with data from the simulation. As 
input, either an extracted feature vector or the filtered wave- 
form is used; as output, binary class labels - bipolar or not - 
are predicted 1311 . The Random Forest and Boosted Trees algo- 
rithms fl51 have achieved the best results for individual sensors 
and clusters of sensors. A Random Forest is a collection of de- 
cision trees. The classification works as follows: The Random 
Forest takes the input feature vector, makes a prediction with 
every tree in the forest, and outputs the class label that received 
the majority of votes. The trees in the forest are trained with 
difl'erent subsets of the original training data. Boosted Trees 



combine the performance of many so-called weak classifiers to 
produce a powerful classification scheme. A weak classifier is 
only required to be better than a random decision. Many of them 
smartly combined, however, result in a strong classifier Deci- 
sion trees are used as weak classifiers in this boosting scheme. 
In contrast to a Random Forest, the decision trees are not neces- 
sarily full-grown trees. For individual sensors, the classification 
error is of the order of 10 % for a well trained model. The com- 
bined results of the individual sensors in a cluster are used as 
new input for training. This method obtains a classification er- 
ror below 2 %. 

4. Conclusion and Outlook 

As shown, the simulation chain is capable of reproducing the 
aspects necessary for acoustic neutrino detection - from the 
generation of the acoustic signal to different detector geome- 
tries and components. Also the deep-sea environment with its 
variable and diverse noise conditions is well replicated. Not yet 
integrated into the simulation is the refraction due to the depth 
dependence of the speed of sound in water Furthermore, the 
whole processing chain for data can be tested as the simulation 
chain employs reconstruction and classification algorithms. The 
algorithms presented are the result of various studies to find the 
most suited ones. Detailed studies about the efficiency of on- 
line filters, effective volume and sensitivity of future large-scale 
detectors are currently pursued. 
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